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Eqnatorially trapped water waves 
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Abstract 

The aim of this paper is to prove that a three dimensional Lagrangian flow which 
deflnes eqnatorially trapped water waves is dynamically possible. This is achieved by 
applying a mixture of analytical and topological methods to prove that the nonlinear 
exact solution to the geophysical governing equations, derived by Constantin in 0 , is 
a global diffeomorphism from the Lagrangian labelling variables to the fluid domain 
beneath the free surface. 


1 Introduction 


In this paper we apply a mixture of analytical and topological methods to establish that a 
recently derived solution dehning Eqnatorially trapped waves is dynamically possible. This 


remarkable solution, derived by Constantin in j5] and given below by equation (2.8), is an 
exact solution of the nonlinear /3—plane governing equations for Equatorial water waves, and 
it is explicit in the Lagrangian framework. The main result of this paper establishes that the 


three-dimensional mapping (2.8) from the Lagrangian labelling domain to the fluid domain 


dehnes a global diffeomorphism— a consequence of which is that the solution (2.8) dehnes 
a fluid motion which is dynamically possible. We achieve this result by hrst establishing 


that (2.8) is locally diffeomorphic and injective, and then we render our results global by 


applying a suitable version of the classical degree-theoretic Invariance of Domain Theorem, 

cf. [ini|22]. 

The solution presented by Constantin in [5] represents a geophysical generalization of 
the celebrated Gerstner’s wave, in the sense that ignoring Coriolis terms in (2.8) recov¬ 
ers Gerstner’s wave solution. The primary importance of Gerstner’s wave is probably the 
fact that it represents the only known explicit and exact solution of the nonlinear periodic 
gravity wave problem with a non-flat free-surface. Gerstner’s wave is a two-dimensional 
wave propagating over a fluid domain of inhnite depth (cf. [31 |1] ES]), and interestingly 
it may be modihed to describe edge-waves propagating over a sloping bed |2l [26]. The 
geophysical solution presented in |5] encompasses Gerstner’s solution, yet it also possesses 
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a number of inherent characteristics which transcends Gerstner’s wave. The solution (2.8) 


is a truly three-dimensional eastward-propagating geophysical wave, and furthermore it is 
Equatorially-trapped— achieving its greatest amplitude at the Equator and exhibiting a 
strong exponential decay in meridional directions away from the Equator. The solution is 
furthermore nonlinear, as is seen from the wave-surface profile, and has a dispersion relation 
that is dependant on the Coriolis parameter. 


Since the solution (2.8) is explicit in the Lagrangian formulation, we may immediately 
discern some qualitative properties of the physical fluid motion. Indeed, an advantage 
of solutions in the Lagrangian framework is that the fluid kinematics may be explicitly 


described [T]. From (2.8) we see that at each fixed latitude the solution prescribes individual 
fluid particles to move clockwise in a vertical plane. Each particle moves in a circle, with the 
diameter of the circles decreasing exponentially with depth. In |5] it was simply shown that 


the solution (2.8) is compatible with the governing equations of the /3—plane approximation 


for Equatorial water waves (2.3)-(2.7). 


The aim of this paper is to rigorously justify that the fluid motion defined by (2.8) is 


dynamically possible. This is achieved by establishing that the solution (2.8) defines a global 


diffeomorphism, thereby ensuring that it is indeed possible to have a three-dimensional mo¬ 
tion of the whole fluid body where all the particles describe circles with a depth- dependant 
radius at fixed latitudes, and furthermore the particles never collide but instead they fill 
out the entire infinite region below the surface wave. In so doing we show that the fluid 
domain as a whole evolves in a manner which is consistent with the full governing equations. 
We note that subsequent to the derivation of Constantin’s solution, a wide range of geo¬ 
physical generalizations and variations to [5] have been produced and analysed, for example 
[6l |71 ISl [HI HH |T71 IlHl lini |2ni El] . It is expected that the rigorous considerations of this 
paper are also applicable to these variants. 


2 The Equatorially trapped wave solution 

2.1 Governing equations 

We consider geophysical waves in the Equatorial region, where we assume that the earth is 
a perfect sphere of radius R = 6378 km. We are in a rotating framework, where the x-axis 
is facing horizontal due east (zonal direction), the y-axis is due north (meridional direction), 
and the z-axis is pointing vertically upwards. The governing equations for geophysical ocean 
waves are given by Euler’s equation with additional terms involving the Coriolis parameter 
which is proportional to the rotation speed of the earth, see 
/ 

ut + uux -k vuy -|- wuz + 2Qw cos <1 — 2^lv sin $ = 

< vt + uvx + vvy + wvz + 2Qu sin <1 = 

wt + uwx + vwy + wwz — 2 VLu cos = 


P 

-\Py ( 2 . 1 ) 

--Pz - g, 
p 
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the mass conservation equation 


Pt + upx + vpy + wpz = 0 
and the equation of incompressibility 

Ux+Vy+Wz= 0 . 


( 2 . 2 ) 


(2.3) 


Here represents the latitude, {u,v,w) is the fluid velocity, H = 73.10“® rad/s is the 
(constant) rotational speed of earth (which is the sum of the rotation of the earth about its 
axis and the rotation around the sun, see m), g = 9.8 m/s ^ is the gravitational constant, 
p is the water density, and P is the pressure. 

We are interested in Equatorial waves, that is, geophysical ocean waves in a region 
which is within 5° latitude of the Equator. Since the latitude is small, we may use the 
approximations sind* ~ <1>, and cos^h « 1, and thus linearising the Coriolis force leads to 


the /3-plane approximation to equations (2.1) given by 


ut + uux + vuy + wuz + 2Ht(; — j3yv 
vt + uvx + vvy + wvz + (3yu 

Wt + UWx + VWy + WWz — 2Hu 

r-ll rvn-lo-l 


= --Px 
p 

= --P 
p 

= --Pz - g, 

p 


y 


(2.4) 


where /? = 2VL/R = 2.28-10 m ^s The relevant boundary conditions are the kinematic 
boundary conditions 

w = pt + upx + vpy on z = ri{x,y,t), (2.5) 

P = Patm on z = g{x,y,t), (2.6) 

where Patm is the (constant) atmospheric pressure, and r/(x, y, t) is the free surface. The 


boundary condition (2.5) states that all the particles in the surface will stay in the surface 


for all time t, and the boundary condition (2.6) decouples the water flow from the motion of 


the air above. We work with an infinitely-deep fluid domain and so we require the velocity 
field to converge rapidly to zero with depth, that is 

(rt, rc) —)• (0,0) as z —>■ —oo. (2.7) 

The governing equations for the /I—plane approximation of geophysical ocean waves are 


given by (2.3)-(2.7). 


2.2 Exact solution 

In this section we present and describe briefly the exact solution of the /3-plane governing 


equations (2.3)-(2.7) which was recently derived by Constantin j5]. This solution describes 


a three- dimensional eastward-propagating geophysical wave which is Equatorially trapped. 
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exhibiting a strong exponential decay in meridional directions away from the Equator, and 
which is periodic in the zonal direction. Equatorially trapped waves propagating eastward 
and symmetric about the Equator are known to exist, and they are regarded as one of 
the key factors in a possible explanation of the El Nino phenomenon (cf. |9l HTj). The 
formulation of the solution employs a Lagrangian viewpoint, describing the evolution in 
time of an individual fluid particle [T]. The Lagrangian positions of the fluid {x,y,z) are 
given in terms of the labelling variables {q, r, s), and time t by 

X = q — sin [k(q — ct)] , 

n 

< y = s, (2.8) 

z = r + pQg ^ 

where k is the wave number, dehned hy k = 2n/L where L is the wavelength, and the wave 
phase speed is determined by the dispersion relation 


c = 


+ fcf/ — n 


and also 


/ «) = 

25 


(2.9) 


( 2 . 10 ) 


determines the decay of fluid particle oscillations in the meridional direction. The labelling 
variables take the values {q, r, s) G M x (— 00 , tq) x M, where rg < 0 is hxed. For every hxed 
s, the system (2.8) describes the flow beneath a surface wave propagating eastwards (in the 


x-direction) at constant speed c determined by (2.9). At hxed latitudes (that is, for s hxed) 


the free surface z = r]{x, y, t) is obtained by setting r = ro(s) in the third equation in (2.8) 
where ro(s) < rg is the unique solution to 


2A:[ro(s)-/(s)] p2krQ 

--rg(s)=" 


2k 


2k 


ro- 


A plot of the free-surface for the wave solution (2.8) is given in Figure below. 

In |5] the author focuses on proving, by explicit computation, that the exact solution 


( |2.8| ) is compatible with the governing equations (2.3)-(2.7). Our aim in this work is to 
prove that it is dynamically possible to have a global motion of the huid domain where, at 
hxed latitudes, the particles move in circular paths with depth-dependant radius. Indeed, we 
prove in our main result Proposition |3.3| that the huid motion dehned by (2.8) is dynamically 
possible, that is, at any instant t, the label map is a global diffeomorphism from the labelling 
variables, {{q,r, s) : g G M, r < rg and s G M}, to the huid domain beneath the free surface 
given by 


{q, s) ~ ’'®’ ’'o(s) + [^(9 “ ct)]^ . (2.11) 

For a hxed latitude s, the surface wave prohle ( 2.11[ ) is a reverse trochoid if rg(s) < 0 and a 
reverse cycloid with a cusp at the wave crest if rg(s) = 0 and s = 0. 
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Figure 1: Exact solution (2.8) 


Fixed s, and given fc > 0 and ro(s) < 0, the curve z = hs(x) given parametrically by 

( 2 . 12 ) 


/ t gro(s)-f(s) pro{s)-f{s) 

I-- 


k 


is a trochoid if ro(s) — /(s) < 0 and a cycloid if ro(s) — f(s) = 0. It represents the curve 
traced by a hxed point at a distance ^ ^ ^ ^ from the center of a circle of radius ^ 
rolling along a straight line without slipping, (see Figure]^. 

Therefore, for a hxed latitude s, the free surface of the huid has the equation z = 
ro(s) + ^ — hs{x — ct) which represents a reverse trochoid propagating to the right with 
velocity c. Since hs is periodic with minimal period 2TT/k then the surface is a periodic wave 
with period 2Tr/k, which concurs with the dehnition of the wave number k. 

3 Main results 


To prove that the motion (2.8) is dynamically possible, it is sufficient to analyse (2.8) for 


the time t = 0, when it takes the form 


X = q 

y = s 


z = r + 


gfc[r-/(s)] 

k 

gfc[r-/(s)] 

~k 


sin (kq ), 
cos {kq). 


(3.1) 
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Figure 2: Trochoid and cycloid curves. 


The case of a general time t in (2.8) is recovered making first the change of variables 


{q, r, s) I—?• (g + ct, r, s), performing (3.11, and finally shifting the horizontal variable x by ct. 


Therefore we can focus on (3.1), and we further note that as q varies by 27rjk, the z value 


reoccurs and x is shifted linearly by 27r//c. Hence, it snfhces to analyse (3.1) on the domain 


= \{q,r,s) : q £ 


27r 


, r < ro and s £ 


In the following result we first prove that the map (3.1) is an injective local diffeomorphism. 


Lemma 3.1. For every fixed t > 0, if rQ < 0 then the map (3.1) is a local diffeomorphism 


from T) = \^{q,r, s) : q £ [O, , r < ro and s £ M} into its image injectively. In the limiting 

case ro = 0, then this result holds for the domain V excluding the cusps at the equatorial 
wave-crests. 


Proof. We remark first that r — f(s) < ro — f{s) < ro < 0, as we see from the definition of 


ro and (2.10). The differential of (|3.1|) at a point {q,r,s) is given by 


1 — e^^'^ cos{kq) sm{kq) —e^^^ sm{kq) 

0 1 0 
—Qk{r-f{s)) g^[Yi{kq) cos{kq) 1 + cos{kq) 


(3.2) 


with determinant 1 — /(^)), As an aside, we note that the time independence of this 


expression implies that the flnid is incompressible and so (2.3) holds, cf. It follows that 


if ro < 0 the Jacobian of (3.1) is non-zero (strictly positive) everywhere, whereas in the 
case ro = 0 the Jacobian is zero precisely at the Equator (s = 0), where the break-down 
in regnlarity corresponds to the appearance of cnsps at the wave-crest as discussed above. 
Therefore, aside from the situation when ro = s = 0, the mapping (3.1) is differentiable. 
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continuous with non-zero derivative, and hence we can apply the Inverse Function Theorem 
to infer that (3.1) is a smooth local diffeomorphism onto its image. 

Let us prove now that (H^) is injective. Let {qi,ri,Si) G P for i = 1,2, and let 
{x{qi,ri, Si),y{qi, ri, Si), z{qi, ri, Si)) be the corresponding fluid particles given by (3.1). First 
of all, if 

(x(gi, ri, si), y(gi, n, si), z(gi, n, si)) = {x{q 2 ,r 2 , S 2 ), y{q 2 ,r 2 , S 2 ), 2 (^ 2 , ^ 2 , S 2 )) 

then Si = S 2 . Thus, we can fix s and then focus on checking injectivity with respect to x 
and 2; in (3.1). Letting ^ = q + ir, then the values of (x, z) in (3.1) correspond to the map 




o-kS{s) 




To prove injectivity, we consider F{^) = ^ -|- 9(0) where ^ = ( 9 , r). Let 0 7 ^ ■^ 2 ) then 
applying the Mean Value Theorem we derive 


|i"(6)-i"(6)l >16-61-19(0)-9(6)1 


> 16 - 61 - max ||T» 9 ^j+(i_s) 5 J| |6 - 61 • 


(3.3) 


Computing Dg in terms of {q,r), yields 

no, . = pHr-m) ( - cos(A: 9) - sin(/c9) \ 
y — sin(A: 9 ) cos{kq) J 


(3.4) 


then ||Il 9 (q,,)|| = /(^)), From (3.3), and considering that 6 = (^i:^*) for * = 1)2, we 

obtain that 

|F(6) - i"(6)l > 16 - 61 - 16 - 61 

= (i-e2^(^'-^(^)))|6-6l) 

where f = max{ri,r 2 }. Therefore, if f — /(s) < 0, T is injective, and we have proved that 
(3.1) is injective. □ 


The following result will be used to prove that (2.8) is in fact a global diffeomorphism, 
cf. [ 22 ], |25|. 

Theorem 3.2. (Invariance of Domain theorem) IfUd is open and F : [7 —)■ M” zs 

a continuous one-to-one mapping, then F : U ^ F{U) is a homeomorphism, and F{dU) = 
dF{U). 


We have already proved in Lemma 3.1 that the exact solution (3.1) gives us a local 
diffeomorphism that is globally injective on F. The result below proves that (2.8) is a global 
diffeomorphism for all t >0, this is, that ( 2 . 8 ) is dynamically possible. 
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Proposition 3.3. For every fixed t >0, if r^ < 0 the map (2.8) is a global diffeomorphism 
from V = {{q,r,s) : g G M, r < tq and s G M} into the fluid domain beneath the free surface 
z = r]{x,y,t). Moreover, if rQ <0 the free surface z = r]{x,y,t) has a smooth profile, and in 
the limiting case rg = 0 the free surface is piecewise smooth with upward cusps at s = 0. 


Proof From Lemma |3.1[ we know that the map (3.1) is an injective local diffeomorphism 
from V = {{q,r, s) : q G [O, ,r < tq and s G M| into its image. To prove that the lo¬ 
cal diffeomorphism is in fact a global diffeomorphism we just have to prove that it is a 
homeomorphism. Indeed, since the hypotheses in the Invariance of Domain Theorem |3.2| 
are satisfied, then the map (3.1) is a homeomorphism. Although it is guaranteed by the 
Invariance Domain Theorem 


3.2 


we can see directly that the map (3.1) sends dT> into 


the boundaries of the image of F. The vertical semiplanes {(0, r, s) : r < vq and s G M} and 
{(27r/A;, r, s) : r < ro and s G M} are transformed by (3.1) in the vertical surfaces 
^Q,y,z)\z < ro(s) -h y G m| and ^{2Tr/k,y,z):z < ro(s) -h y G m| 

respectively, and the horizontal semiplane {((7,ro,s) : 0 < g < 2TT/k and s G M} becomes 
part of the reverse trochoid if ro(s) — f(s) < 0, which is smooth, and it becomes part 
of the reverse cycloid if vq^s) = 0 and s = 0, which is piecewise smooth with upward cusps. 


We have proved that (3.1) is a global diffeomorphism map from F into its image if 


To < 0, with singularities occurring when tq = 0 and s = 0. Since the full system (2.8) can 


be recovered from (3.1) by making the change of variables {q, r, s) i—)• {q + ct, r, s), and finally 


shifting the horizontal variable x by ct, it follows that (2.8) is a global diffeomorphism from 
V = {(y, r, s) : g G M, r < ro and s G M} into the fluid domain below the free surface. □ 
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